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TECHNICAL NOTE 39 k 1 


COMPARISON OF CALCULATED AND EXPERIMENTAL LOAD 
DISTRIBUTIONS ON THIN WINGS AT HIGH 
SUBSONIC AND SONIC SPEEDS 
By John L. -Crigler 

SUMMARY 


A method for calculating the aerodynamic loading on a wing in combi- 
nation with a body is presented. The method is similar to that used by 
Falkner for wings in incompressible flow, except that the present method 
directly relates the downwash velocity to the lift on the wing in a com- 
pressible medium for any Mach number. 

Calculated results are compared with experimental data for two wing- 
body configurations throughout a range of Mach number up to 1.0. Both 
wings were swept back 45° with an aspect ratio of k . One of the wings 
hafl a taper ratio of 0.6 and NACA 65 AOO 6 airfoil sections and was tested 
with and without twist. The other wing had a taper ratio of 0.15 and 
NACA 6kA206, a = 0 airfoil section at the fuselage boundary, fairing to 
an NACA 6kA203, a = 0.8 (modified) airfoil section at the midspan. The 
sections of the wing from the midspan to the tip were NACA 6kA203 , a = 0.8 
(modified ) . 

The magnitude and the distribution of spanwise loading of the cal- 
culated data are in good agreement with experiment up to a Mach number 
of 0 . 95 , and for the highly tapered wing the agreement of the calculated 
spanwise load distribution with the experimental distribution is still 
good up to a Mach number of 0.98* 


INTRODUCTION 


Tn recent years the calculation of span loading on wings for wing- 
body combinations has become of increasing importance for airplane 
designers. The great number of airplanes with widely different types of 
plan form, designed to fly at high Mach numbers, makes it virtually 
impossible to obtain experimental data on the aerodynamic characteristics 
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over the complete flight range for all wings of interest. Thus, there 
has been considerable interest in calculation methods, particularly in 
the calculation of the distribution, or shape, of the span loading. 

The purpose of this report is to show to what degree span load 
distributions as calculated by linear theory agree with experimentally 
measured span load distributions, particularly in the speed ranges 
approaching sonic conditions. In making this study, calculated span 
load distributions are compared with experimentally measured values for 
some thin wings for wing-body combinations throughout a Mach number range 
up to M Q = 1.0. 

The method used in making the calculations for the present study is 
similar to that used by Falkner in reference 1 for treating wings in 
incompressible flow and to that used in reference 2 for treating oscil- 
lating finite wings in subsonic compressible flow. The present method 
directly relates the downwash velocity to the lift on the wing in a 
compressible medium for any Mach number up to M = 1.0. As in Falkner' s 
method for incompressible flow, the loading is assumed to be expressible 
as the sum of terms involving certain unknown coefficients, appearing 
in the lift distribution. By assuming the downwash to be known at an 
appropriate number of points, a system of linear algebraic equations is 
obtained from which the unknown coefficients can be determined. 

In order to calculate the load distributions of the wing in combina- 
tion with a body use was made of the method suggested in reference 3, 
wherein the boundary conditions on the fuselage are assumed to be ade- 
quately satisfied by a vortex image system. 


SYMBOLS 

A aspect ratio, ks 2 /S 

I 

a^^ coefficients in expression for lift 

C L lift coefficient, L/qS 

c section chord, ft 

c r root chord, ft 

S 

c &v . average chord, ft 
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c t 

c l 

c m 

D 

d 

E 

F, f 

G, 6 
%i,v 
K 

L 

l 

Al 

M 

Mo 

m 

m' 


tip chord, ft 

section lift coefficient 

section moment coefficient 

fuselage diameter, in. 

arbitrary distance, ft 

total downwash factor 

downwash factors for trailing vortex 

downwash factors for bound vortex 

chordwise replacement load 

kernel function 


lift 


, Jl dy. 


lb 


section lift 


, J' Al dx. 


lb 


lift at any point, lb 
moment, ft-lb 
free-stream Mach number 

tangent of angle of sweep of local chord line 
tangent of angle of sweep of midchord line 


qL dynamic pressure, V^, lb/sq. ft 

r radius of an infinite cylinder representing body, ft 

r dimensionless radius with respect to 

S wing area, sq. ft 

s wing semispan, ft 

a^. semi span of vortex 
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V 

w 

X 


x,y,z 

x,y 





air stream velocity, ft/sec 
induced vertical velocity, ft/sec 

distance of lift point from midchord line along X-axis, 

Xm " x b 

longitudinal, lateral, and vertical coordinates, respectively 

x y 

dimensionless coordinates, ^ 

coordinates of points along -midchord line 

wing section center of pressure 

longit udina l center of pressure for wing-body configuration 


a 


angle of attack, deg 



M 


5 


y 

T) dimensionless coordinate, j 

9 = COS_1 Tpk 

X taper ratio, c^./c r 

A sweep angle of quarter -chord line, positive for sweepback, 

deg 

p density of air, slugs/cu ft 


Subscripts: 

a index of control points 

b index of lift point or vortex 

Superscripts : 

(i) downwash factor of image vortex 

' downwash factor of vortex located at point (x^ , -y^J 



NACA TN 39^1 


5 


METHOD 


The method used in making the calculations is as follows: 

Figure 1 shows a wing diagram and the coordinate system employed. 
In reference 4 the boundary -value problem for the lift distribution on 
a finite wing is reduced to an integral equation of the form: 


w (Wa) = & 



x b^b - ^b 


(i) 


In this expression w^x a ,y a ') represents the displacement velocity at 
any point ^x & ,y a ^ , or the local downwash velocity at point ^x a ,y a j, and 

may be assumed to be known. The value Al represents the unknown lift 
distribution and the kernel K represents the known vertical velocity, 
or downwash, at point (x a ,y a ^ associated with a pressure doublet of unit 

strength located at any other point ^x b ,y b ). In the case of incompressible 
flow, K is found to be a function only of the distances (x a - x b ) n.wfi 
(y-b " ya)- 111 tlie case of compressible flow, K is a function not only 
of the distances (x a - x b ) and (y b - y a ) but also of the Mach number M Q . 

The integral equation for this case may be written as (Bee eq. (15), 
ref. 4) 




(V *aY 


1 rf* 


x a- x b 


J( x a- x b) 2 + p 2 (yb-y a ) : 


|dx b ^b 
( 2 ) 


where S denotes the surface of the wing projected on the xy -plane. 

The problem now is to determine Al ( x b ,y b ) so as to satisfy equation (2) . 

For arbitrary wings, it does not appear possible to obtain exact solutions 
to this integral equation; therefore, approximate methods must be resorted 
to. The method employed is to assume that the lift Al (x b ,y b ) can be 

expressed in the form of a series of loading terms having both spanwise 
and chordwise variables, as 
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8 P v 2 

Al (V^b) = — tan a 


nP 7 


*b 



+ 


sin 





+ 


sin 


29 





( 3 ) 


where 


cos 0 


x n ~ x b _ X 
c/2 ” TJI. 


(o < 0 < it) 


as shown in figure 1. The value 0=0 corresponds to points on the wing 
leading edge and 0 = it corresponds to points on the wing trailing edge. 

In equation (3), s denotes the semispan of the wing so that this 
expression for the lift denotes an elliptic type of loading with modifi- 
cations given in terms of the variable y^. The first term of the dis- 
tribution, cot implies a -y- type of singularity aft the wing leading 

edge, and each term of the series goes continuously to zero at the wing 
trailing edge. Substituting the expression for the lift given in equa- 
tion (3) into the integral equation (2) expresses the total local down- 
wash w(x a ,y a ) as a sum "terms involving the coefficients a^^ as 

the unknowns with all other quantities in the equation known. For 
example , 




il 1 

j i 

h 

i 

£ 

1 

; J (y-b - y a ) 2 

sj ( x a - x b) 2 + e 2 (yb - y a ) 2 


sin 8 d0 
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or 



The problem now is to evaluate the integrals in equation (4) in order to 
determine the coefficients It should be noted, however, that 

equation (4) consists of a series of elliptic integrals which cannot be 
evaluated in closed form. The integrals can nevertheless be accurately- 
performed by approximate or numerical procedures for any preassigned 
values of ( x a ,y a )- Since the total downwash velocity w ( x a ,y a ) I s 

assumed to be known, it may be seen that the integrations yield equations 
for 'w(x a ,y a ) in terms of the unknown coefficients appearing in 

the lift distribution. Thus, by assuming the downwash velocity to be 
known at an appropriate number of points, called "control points," a 
system of linear algebraic equations is obtained from which the coeffi- 
cients an jm can be determined. It is significant to note that, when 
a is determined, the lift is expressed in a perfectly continuous form 
except for the allowable singularity along the leading edge. 


The Spanwise Integration of the Downwash 

In evaluating the integrals in equation (4), either the chordwise 
or the spanwise Integration may be performed first. The procedure herein 
begins with the spanwise Integration) that Is, with an evaluation of the 
integral 



8 


me A TN 3941 


ff(Wa) " 6 ^ AZ(9,y b )K(x a ,y a , M 0 , 9, y b ) dy b (5) 

where x a , y a , M 0 , and 8 are assumed to he fixed. In order to evaluate 
the integral along a curve of constant 0, let this curve be divided 
into small arc lengths Ad by dividing the span length 2s into small 
elements of span, of width 2s y , and let x^,y^ represent the midpoint 

of the pth span element as indicated in the following sketch. 



The integral, equation (5), may now be written as a sum of integrals, 

s/s Y r y^+s v 

^0'( x a->^a) * n^T sin 9 zH / ^( e .»yb) K ( x a>ya> Moi 9 > ^^b 

p=l J y^-s v 


( 6 ) 

The approximation is made on the assumption that the element of span of 
width 2s v can be taken small enough, compared with the semi span s, 
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that the value of the loading A 1 , y^) at the midpoint of the element 

adequately represents the value at any point of the element. This span- 
vise approximation may be carried out to any desired degree of accuracy 
so that in the limit the feature of a continuous load distribution is 
retained. 

Thus, the value of — may be approximately -written as 


s/s* v* 

^( x a.y a ) w rj 2Z 2s v Al ( 9 ; yp)E(3^, M Q , 0, y£-y^) sin 6 (7) 

^ (i=l 


where 


_ o^ +s v 

E ( x a* ^ 0 ’ yp-^a) = / K ( x a^a^ M 0 , 9, 

J yp-s v 


yb) 3 ^ 


(8) 


The integration required to obtain the function E in explicit form 
can readily be carried out with the results expressed in terms of the 
distances (x a - x^) and (y^ - y a ) and the free-stream Mach number Mq. 

(See, for example, ref. 5 where the expression for E is derived by using 
a different notation. ) The value of E is identical with the downwash 
function associated with a horseshoe vortex where the bound leg is of 
length 2sy with its midpoint located at x (i ,y 1 _ l . As will be discussed 

later, it is necessary to apply a correction to the downwash due to the 
bound leg of the horseshoe vortex to account for the increase in vertical 
velocity caused by the fuselage. For this reason the contributions 
associated with the bound legs and trailing legs are employed separately. 
The derivation of E in this separated form is given in appendix A. 


The Chordwise Integration of the Downwash 

After having performed the spsnwise integration, the next step is 
to effect the chordwise integration j that is. 


p s/s v nit 

W ( x a^a) = 77 2s v 2Z J Q Al ( Q > yp) E (^ M o> 0 ; V*y^)sin 9 d6 


(9) 
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It must be remembered that the lift distribution given by Al ^0 , y^j 

in equation (9) is an unknown quantity, involving certain coefficients 
a^ m which have to be determined. The integrations in equation ( 9 ) naay 

be effected by evaluating the integrand for fixed values of x a , y a , M 0 , 
and y at small intervals of 0 and then making use of some curve- 
fitting scheme of integration such as Simpson’s rule. It should be noted, 
however, that since E is equivalent to the downwash function associated 
with a horseshoe vortex, equation ( 9 ) is the basic equation of Falkner’s 
procedure for treating the lifting surface. Falkner, in lieu of using 
fima. l 1 spacings in the chordwise direction, devised a scheme to approximate 
the integration that involves evaluating the integrand in equation ( 9 ) for 
only a few discrete values of 0. Falkner used the method for the case 
of a steady wing in incompressible flow, but the method can easily be 
extended to steady wings in compressible flow and to oscillating wings 
in either compressible or incompressible flow. (See, for example, ref. 2.) 
The method of chordwise integration, as devised by Falkner, consists of 
expressing Al(Q , y^) in terms of Bq m for only a few discrete values 
of 0, denoted herein by 0 V (v = 1, 2, 3, • • • Q) and employing weighting 
factors G^v associated with each value of v. Although there is no 
apparent justification for using this method of chordwise integration for 
Mq = 1, calculations have been made for this condition herein. 

The method of determining the weighting factors v is discussed 
in the next section. 

It Is convenient to locate the control points so that they fall along 
fixed percentage span and chord positions; that is, along constant tj _ and 
constant 0 curves. The spanwise locations of control points are chosen 
so that they lie midway between the legs of the horseshoe vortex. The 
chordwise locations are chosen so that the control points fall midway 
between the lines of constant 0, which means there must be at least one 
more discrete value of 0 than chordwise control points. Since one equa- 
tion is formed at each control point there, should be the same number of 
control points as coefficients retained in the equation for lift. 

In the present treatment, the number of both spanwise and chordwise control 
points are made to equal the number of spanwise and chordwise terms retained 
In' the series, nine. 

Now, at any control point with coordinates x a ,y a , let G n ^ v repre- 
sent the weighting factor for the nth term of the spanwise loading distri- 
bution and the vth value of 0. The ^ith integral in equation ( 9 ) may then 
be written 


w (Wa) 



Al (8, y^) e(x bj Mo; 8, y^-ya^sln 8 40 


« i. 

kit 



o + a nA n 4- . ._. + 1^, e vJ yjr-yj)- (lo) 
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Determination of the Weighting Factors 

Hie weighting factors G n ^ v are determined with the use of theory 

for a two-dimensional wing by setting up a system of algebraic equations 
for the lift ahd the downwash involving these factors. Each chordwise 
term in the lift series requires a system of equations j but each set of 
equations is obtained in the same manner. The details are therefore 
presented only for the n = 0 term; that is, for the term containing 

cot %. 

The first equation to be considered equates the integral of the 
continuous loading to the sum of the weighting factors; namely, 



= G 0,l + G 0,2 + * * * + G ( 


’•« = / 0 


cot sin 0 d9 = | (ll) 


Each of the other (Q-l) equations in the system involves one chord- 
wise location of control point (midway between lines of constant 0) and »n 
the discrete values of 0. These equations are derived by the use of the 
two-dimensional kernel or downwash function. The downwash at the control 
point involved, as produced by the sum of the loads on each discrete value 
of 0 is equated to the downwash produced by the continuous loading in 


two-dimensional flow; that is, if x^, 


: (Q-1) 


denotes the 


(Q - l) chordwise control-point location (x is midway between lines of 
constant 0), then for x = x-. 


L" O 0 , v K(x lJ M OJ 0y) - 0 0)1 K( Xl , Mo, 6 x ) + M„, 0 2 ) 


0o] + • . . + 


^ 6 m) 


1# 


= / K^X]_, Mq, 0)cot ^ sin 0 


^0 
= jtp 


d0 


( 12 ) 


These equations are independent of Mach number and, therefore, are iden- 
tical to the equations for the incompressible case. The weighting factors 
for various values of 0 y have been tabulated. (For example, see ref. 1. ) 
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The equations for the first term of the lift series ^ the term con- 
taining cot for four values of 0 corresponding to the 12. 5- j 37*5-; 
62. 5-j and 87.5-pencent-chord lines are 


G 0,l + G 0,2 + °0,3 + G 0,4 = | 

8G 0,1 “ 8G 0,2 " 5 G 0,3 “ 5 G 0,4 = * 
f G 0,l + 8G 0,2 ■ 8 g o,3 " 3 G 0,4 = 
f G o,i + 3 g o ,2 + 8 g o,3 " 8 g o,4 = 


(13) 


Similarly, four equations to determine & n>v for any other value of n 
are obtained by replacing the load term containing cot ^ in equation (ll) 


by the term containing sin n0. The weighting factors are tabulated in 
table I for four values of 0 corresponding to the 12. 5-> 37* 5~> 62.5-, 
and 87. 5-percent -chord lines. 


Determination of the Downwash 

After G n v is determined for a given choice of 0 -values, the 

expression for the total downwash at any particular control point is 
given with the use of equations (9) and (10) as 


V 


(xa^a) 


1 



%ii + 


7* 

+ a n,*0 E (^' V *v> 


Since the total downwash is known at the chosen control points, equa- 
tion (l4) leads to a system of linear algebraic equations from which the 
coefficient can be determined. 


S’ 
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An example showing the formation of the downvash equations for the 
84 -vortex pattern with four 8 -lines is given in appendix B. 


Effect of Presence of Fuselage 

In the analytical treatment of the wing-fuselage combination, the 
fuselage is represented by an infinite cylinder, which is replaced by 
the images of the wing vortices placed inside the fuselage as suggested 
in reference 3. As indicated in reference 3, the boundary conditions 
are satisfied exactly only infinitely far behind the wing. The vortices 
used to represent the fuselage induce a downwash on the wing which must 
be taken into account. This additional downwash may be approximated by 
a correction to the downwash due to the bound vortices as discussed in 
reference 3. The corrected downwash is then given as 


Corrected downwash = Bownwash [1 + 


MODELS 


Calculations have been made for two wing-body configurations. One 
wing, designated wing A, was swept back 45°, with an aspect ratio of 4, 
a taper ratio of 0.6, and KACA 65 AOO 6 airfoil sections parallel to the 
plane of symmetry. Details of wing A in combination with the body are 
given in figure 2. Experimental data used in the comparisons for wing A 
were obtained from tests made in the Langley 8 -foot transonic tunnel. 

Both experimental and calculated results were obtained for this Wing 
without twist and for the wing with a spanwise twist variation from the 
fuselage boundary to the wing tip as shown in figure 3- ThiB spanwise 
variation was produced by twisting the wing sections about the wing 
quarter-chord line in planes parallel to the vertical plane of symmetry. 

In making the calculations the body in all cases was treated as an infi- 
nite cylinder, as indicated by the dashed lines in figure 2 . 

The other wing, designated wing B, was swept back 45°, with an aspect 
ratio of 4, a taper ratio of 0.15, and UACA 64A206, a = 0 airfoil sections 
at the wing midspan, fairing into KACA 64A203, a = 0.8 (modified) airfoil 
sections at the 0.5 semispan and retaining these sections from the midspan 
to the wing tip. Details of wing B in combination with the body are given 
in figure 4. These calculated data are compared with unpublished experi- 
mental data obtained in the Langley 8 -foot transonic tunnel. 
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COMPARISON OF EXPERIMENTAL AND CALCULATED RESULTS 


Figure 5 shows a comparison of the calculated spanwise loading of 
wing A alone with that of the same wing in combination with a body. The 
details of the wing-body configuration are shown in figure 2. The cal- 
culations are for an angle of attack of 4-.0°. Figure 5(a) shows the 
results for incompressible flow (Mq = 0) and figure 5(h) shows the results 
for compressible flow for Mq = 0.9« The only significant effect of the 
body is to reduce the loading on the wing adjacent to the body and for 
these two cases the effect was less at Mq = 0.9 than for the incompress- 
ible flow case. The loadings on the more important outboard sections are 
not significantly affected. In all remaining calculations the effect of 
the body is included for comparison with experimentally measured data 
since no experimental data have been obtained with the wing alone. For 
these and all subsequent calculations, the wing is assumed to have a 
fixed geometry; that is, no deflection or bending under load. For the 
wings tested in this report, deflections under load were very small at 
a = 4-.0°, so that any discrepancies introduced in the calculations due 
to deflection or bending were small. However, in many cases, this effect 
is important, and if so, must be included in the calculations. 

Figure 6 shows a comparison of the calculated and experimentally 
measured spanwise loading plotted against the wing semispan for Mach 
numbers of 0.6, 0.95; and 1.0 for wing A in combination with a body. 

The data are for an angle of attack of 4.0° and for the wing without 
twist. The spanwise distribution of lift, on which the wing spanwise 
shears and bending moments depend, are generally in good agreement, both 
in shape and magnitude. It may be noted that, for Mq = 1.0, there is 
a drop in the measured lift adjacent to the body and also an increase 
near the tip which results in a change in the shape of the loading curve; 
that is, an outward shift in the lateral center of pressure which is not 
found in the calculated results. 

Figure 7 shows a comparison of results for the same wing-body con- 
figuration and the same Mach number range as shown in figure 6; however, 
in this case the wing sections were twisted about the quarter-chord line 
to produce the spanwise twist variation shown in figure 3. These data 
are also for an angle of attack of 4.0°, the angle of attack being measured 
from the body center line. The agreement of theory and experiment, here 
again, is good and the same trends with Mach number as in figure 6 are 
noted. It thus appears that the changes in span loading which result 
from twisting due to aeroelastic bending or wing deflections can be pre- 
dicted at these Mach numbers. 

Figure 8 shows calculated and experimental spanwise loading coeffi- 
cients for wing B. The details of the wing-body configuration are shown 
in figure 4. The results are for an angle of attack of 4.0° and are shown 
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for Mach numbers of 0.80, 0 .94, 0. 98, and 1.0. The agreement between 
theory and experiment for Mach numbers of 0.80 and 0.94 shown in fig- 
ure. 8(a) and for Mq = 0.98 shown in figure 8(b) is very good. This 
wing was not tested at a Mach number of 1.0, the highest Mach number for 
which calculations have been made, but experimentally measured data were 
obtained at Mo = 1.03. These experimental data fell below the curve 
for Mq = 0.98 and there was no appreciable difference between the 
experimental data for Mq = 0.94 and Mo = 0.98. The shape or distri- 
bution of the loading for the calculated data at Mq = 1.0 was essen- 
tially in agreement with the distribution of the experimental curves for 
Mq = O .98 and Mq = 1.03. Thus it appears that for this wing there is a 
an appreciable increase in the magnitude of the loading on going from 
Mq = O.98 to Mq = 1.0 for the calculated data that would not be found 
experimentally. 

For the cases that have been calculated it may be concluded that 
the magnitude and the distribution of spanwise loading calculated for 
the thin wings under consideration are in good agreement with experiment 
up to a Mach number of 0.95, and for the highly tapered wing (wing B) the 
agreement is still good up to a Mach number of O.98. 

Figure 9 shows a comparison of the calculated and experimentally 
measured chordwise center -of -pressure location for several spanwise sec- 
tions for wing A over the range of Mach number tested and figure 10 shows 
similar results for wing B. The experimental data in figure 9(a) for the 
untwisted wing for Mq = 0.8 and M Q = 0,9 were taken from reference 6. 
On both wings the agreement over the inboard section is good up to a Mach 
number of 1.0. However, there is a sharp rearward shift in the chordwise 
center of pressure location shown at the 0.8 and 0.95 semispan stations 
for the experimental data occurring between Mq = 0.95 and M 0 = 1.0 
that is not indicated by the calculated data. With regard to the chord- 
wise pressure distribution good correlation would not be expected, par- 
ticularly over the outboard sections where there might be a tendency for 
the flow to wash out towards the tip. 

Figure 11 (a) shows a comparison of the calculated and experimental 
longitudinal center-of -pressure location relative to the leading edge of 
the root chord versus Mach number for wing A with body without twist. 

The calculated data are for a constant angle of attack of 4.0°. The 
calculated lift coefficients, varying from = 0.260 at Mq = 0.8 to 
Cl = 0.306 at Mo = 1.0 are shown on the figure. The experimental data 
are for lift coefficients of 0.2 and 0.4. The calculated longitudinal 
center of pressure for the wing -body combination is approximately at the 
trailing edge of the root chord for this wing without twist and is seen 
to shift slightly rearward with an increase in Mach number. The shift 
in the longitudinal center of pressure for the experimental data is 
somewhat greater as Mq approaches 1.0. This greater shift in the 
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longitudinal center of pressure for the complete wing would be expected 
from the sharp rearward shift in the chordwise center-of-pressure loca- 
tion for the outboard sections shown in figure 9* 

Figure 11 (b) shows a comparison similar to figure 11 (a) for wing A 
with body and with twist in the wing. Again the calculated data are for 
a constant angle of attack of 4.0° and the calculated lift coefficients 
for several Mach numbers are shown. In general, the same trends are 
shown, except in this case the longitudinal center-of-pressure location 
for the experimental data at = 0.1 was about 7 percent of the root 
chord behind that for a of 0.2. Other experimental data for this 

configuration, not presented herein, showed no appreciable shift in the 
center-of-pressure locations for lift coefficients of 0.6 and 0.8 from 
that at Cl = 0.1. 

Figure 12 shows a comparison of the calculated and experimental 
longitudinal center-of-pressure location relative to the leading edge 
of the root chord against Mach number for wing B. The calculated data 
are for a constant angle of attack of 1.0°. The experimental data are 
only given for a lift coefficient of 0.1 since the calculated lift coef- 
ficients were near 0.1. The same trends as noted in figure 11 are shown 
for this wing -body combination. 


CONCLUDING REMARKS 


A detailed method for calculating the span load distribution of a 
wing in combination with a body is presented. A comparison of the _ 
calculated load distributions with experimental results indicated that, 
for the limited number of cases calculated, the magnitude and distri- 
bution of spanwise loading calculated for these thin wings are in good 
agreement with experiment up to a Mach number of 0.95 and that the span- 
wise loading calculated for the highly tapered wing is in good agreement 
with experiment up to a Mach number of 0.98. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 11, 1958. 
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APPENDIX A 

DERIVATION OP DOWNWASH FACTOR E 


It was pointed out in the text of the report that the method of 
effecting the spanwise integration of loading employed led to a system 
of horseshoe vortices. From considerations in reference 3, the wing- 
body combination is represented by a set of vortices located on the wing 
with image vortices located inside the body. Since a correction iB 
applied to the downwash associated with the leg of the vortex, as out- 
lined in reference 3, it is desirable to calculate the contribution due 
to the trailing legs and bound leg separately. The downwash in this 
separated form is calculated with the use of the Biot-Savart law. 

The contribution of the trailing IjLegs of the vortex at (x^^y^) and 
of width 2s v to the downwash factor at point (x a ,y a \ is given by 


1 

( x a " x b) + \l 

( x a - x b) 2 + P 2 (yb - y a + s v) 2 

y b " y a + s v 
1 

/( x a - X 1 
( x a - x b) + J 

3 ) 2 + p2(y b . y a * S t )2 

—1 

( x a - x b) 2 + P 2 

CVI 

to 

1 

s ^ 

t> 

CO 

1 

£ 

1 

£ 

\fl*e. - x b) 2 + P 2 (^b * Ta - s v) 2 J 


(Al) 


or referring all distances to the vortex semispan s v 


1 1 

(% - x b) +\l ( x a - x b) 2 + P 2 (yb - y a + i) 2 

- y a + 1 s v 

\i{\ - ^t) a + ^ - ^ +i ) 2 

i i 

( x a - x b) + v( x a - x b) 2 + P 2 (y-b - - !) 2 

*b ‘ ^a : 1 s v 

_ /( x a - x bf + P 2 (^b - y a - 1 ) 2 


(A2) 


The contribution of the bound leg of. the vortex at (x^y^ and of width 
2s v to the downwash factor at point Xa>y a which is increased by the 


factor 


(“$ 


(see ref. 2) is given by 
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The contribution to the downwash factor at point (x a ,y a ) due to the vortex 
at (*^-7*,) f° r a wing with symmetrical loading, Al = A£ (x^, 

is given by 




Each discrete vortex located on the ving is considered to have an 
image vortex of the same strength inside the fuselage that must be con- 
sidered in summing up the dowrrwash at any control point. A horseshoe 


vortex of -tfidth 


located, outside of a cylinder of radius 


an image horseshoe vortex of width (- —I within the cylinder. Then 

Yl y 2 j 

the contribution to the downwash factor at point (x a ,y a ) due to the image 
inside the fuselage of the vortex at x- b ,y- b is given by 
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The contribution to the dovnwash factor at point Xa,>y a due to the image 
inside the fuselage of the vortex at , -y^ is given by 




20 


NACA TN 391+1 



The total dovnwash factor at point x a ,y a due to the vortex at 
x^y^, its companion vortex at ^,-y-^ (on a symmetrical wing), and 
their images inside the fuselage are given by 
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APPENDIX B 

APPROXIMATE INTEGRATION IN THE CHORDWISE DIRECTION AND 
FORMATION OF THE DOWNWASH EQUATIONS 


In the formation of the downwash equation, there must be one control 
point for each equation. The location of the control points on the wing 
depends on the layout used, and the number and optimum locations of the 
points must be determined by experience. Falkner found the 84-vortex 
pattern with nine control 1 points located as shown in f igure 13 to be 
satisfactory in computing the load distribution of a s ymme trical swept- 
back wing in incompressible flow with the 328-vortex pattern giving only 
slightly better accuracy. Therefore, the 84- vortex pattern, or the 
equivalent, for the wing-body combination which divides the wing into 
intervals of 0.1 semispan and four chordwise load lines has been used 
for all calculations in this report. It should be noted, however, that 
it is quite probable that placing the control points nearer the trailing 
edge would give better results, particularly as M Q — >1.0. Both the 
spanwise and chordwiBe distribution of control points are the same as 
the spanwise and chordwise terms retained in the equations. 

The equations for determining are formed as follows: 

The coefficient for ^ corresponding to the cot ^ term for 

the 84-vortex solution (control point 1, fig. 11) is given by (see 
table I for values of G n ^ v ) 


0.2734 22 E-l/8 n/i - n 2 + 0.1172 22 e 3/8\A - * 


0.0703 22 E 5 / 8N /l - tj 2 + 0.0391 2Z E 7/8N I 1 - T l 2 


where the subscripts l/8, 3/8, 5/8, and 7/8 on the downwash factors refer 
to the vortices on the 12. 5-, 37-5-, 62.5-, and 87. 5-percent -chord lines, 
respectively. 
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The coefficient for q corresponding to the sin 0 term is 
given by 


o.o 488 £2 e ]_/8 v/ 1 - n 2 + 0.0762 ^ E^g ^1 - tj S + 

0.0762 ^ V8 JI - t) 2 + o.o 488 E7/ 8 yi - n 2 

The coefficient for a 2 ^ 0 corresponding to the sin 20 term is 
given by 


0.0732 \ fa -~ J 1 “ 1 2 + 0,038:L H e 3/8\A " r ‘ £ ' 

0.0381 2Z E 5/a n/ 1 ' ^ - °-° 732 X[ =7/8^ ‘ ^ 

For more than three chordwise control points it would be necessary to 
include terms of sin n0 where n — >2. 

Additional spanwise terms in the equation are determined by adding 
further spanwise loading functions . The coefficient for aQ 2 corre- 

0 

sponding to the second cot ^ term is given by 

0.2734 E 1 /8 t > 2 n/i " T » 2 + °* 11T2 51 ^/s^ 2 'i 1 “ 4 2 + 

0.0703 ^ E 5 / 8 ti 2 v/l - t} 2 + 0.0391 E 7/8 n2 'J 1 ~ 7)2 

Similarly, the coefficients for the other unknowns corresponding to the 

cot the sin 0, and the sin 20 terms are found. For each equation 

formed the downwash velocity is equated to the local slope of the wing 
section; that is, w/v is equated to tan a or, for the wing divided 
into intervals of 0.1 semispan, each equation is set equal to 0.025 to 
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allow for the ratio — = kO. In a similar mann er, equations are formed 

s v 

for the other control points, and, for the case of the wing without twist 
or camber, all equations for all control points are set equal to 0.025. 

If twist or camber is incorporated in the wing, the constant for each 
equation will depend on the twist or camber at the particular control 
point j that is, the local slope of the wing section. 

For the wing -fuselage combination the effective angle of attack 
must be corrected due to the vertical velocity component induced by 
the fuselage. If it is assumed that the downwash induced by the bound 

*p»2 

vortices in the plane of the wing is increased by the factor 1 + 

y a 

as in reference 3 , the effective angle of attack of the wing is equal to 
the geometric angle plus the angle induced by the fuselage 


or 
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APPENDIX C 


AERODYNAMIC CHARACTERISTICS 


The lift (see eq.. 3) for any point is given by 


SpV 2 


Al = s tan 


ajrr 


cot l( a o,o + ^0,1 + A 


,2 * * •) 


sin 0 + T l a i #1 + *1 ®i # 2 * * *) + 

sin 20 + n« 2 A + n% >2 • • •) + • • • 


(Cl) 


j . agi „ tan f /s L |(a 0> 


or the total lift at any chord is 

'c/2 

-=/2 l 

sla e(a lj0 + T)ai,l + 1^1,2 +•••) + 

aln 29 ( a 2,0 + ^2,1 + lS,2 + ' ‘ + ' ' j 


0 + + a 0,2 + 


•) 


ax 


(C2) 


but 


and 


Then 


X = ^ cos 9 


ax = - i sin e ae 


/ c / 2 pC 

cot % dX = - / cot & sin 9 d0 = & 

c/2 
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f c /2 

/ sin 6 dX = f 
J -c/2 k 


and 


pc/2 

-c/2 


sin 20 dX = 0 


Therefore, equation (C2) for a symmetrical wing becomes 


l = 8pV*B tan c^l -\ 2 [§f 0 ,2 + *1%^ + + 

t ( a X,0 + 1 2a l,2 + 


or the section lift coefficient is 


c z = 


SpV^an a 
JpV 2 c 


sf ~- ’l 4 2 ( a 0,0 + ’>**0,2 * + 


(*1,0 + 1 ' 2a l,2 + iS.lt) 


but the total lift on the wing is 


■'S 


L = f l dy 
J -s 


or 


c L - 


s / i dt) 
'-1 

|pV 2 S 


(C3) 


(Ch-) 
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and 



Cl = £b__ q tan a I Jl - t\ c 
S ^-1 


£ 

4 


^/a + Tl 2 a^ ^ + T^a 

2 l 0,0 0,2 


dr} 


0,4) 


( a !,0 + *\,2 + * \,| 

l(l a 0,0 + f a i,o) + 8(2 a o,2 + IT a i,2^ 

JC fa It \ 

16\2 a o ,4 4 \ k ) 


Cl =' 4A tan a 


where 


4s 2 


S 


Then 


C L = ^ - (l6e. + 8a + 4a + 2a + 2a . + a 

l 6 \ 0,0 1,0 0,2 1,2 1> 


The slope of 
da 


the lift curve Is 

^-^.6a 0 ,o + 8 ^l,o + ^ a 0,2 + 2a l,2 + 2a 0,4 + a l,4^ 


The ratio of the section lift coefficient c^ at any section 
total lift coefficient Cj, is 


Cl _ l6Sy/l - 
Cl scst 

fa(«0.0 + l a °0.2 * l 4 °0.it) + (»i.o •*■ Al.2 - 1 - 

[ l6a 0,0 + **0,2 * Z \k + \o + **1,2 + a l,^ 


> (05) 


^ (06) 


(07) 

to the 


(C8) 
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The section moment coefficient about the midchord is given by 


c m = 


W tan a, 


Z 1 “ ^ ^ i - (i - |a 0 ^ 0 + 


+ ^ a A k + 


2^2, 


0 + 71 a 2 ; 2 + * a 2 


and the center of pressure of any section is given by 
c 

C P _ m 
c c z 

_ c av [l- (1 -^[iao.o + ^ao.a* l‘ ta 0.1t)' l '?( a g.Q' l ' , 1 gB 2.g' l 'l‘ ta 2,l|)l 

[ 2 ( a 0,0 +, ' 2a 0,2 + ’' 1,a 0 ) l) t ( a l,0 + ’lS,2 + "\,lJ 


The moment for the complete wing about a line through the leading 
edge (x = 0) of the center chord (minor axis) is given by 


_ , 8pV 2 B 2 (?T + £ _ ~2 M + ^a o>2 + 


n a n i , + 


X3,k 2 


*1,0 , 2 a l,2 , 4 a l,4 


dT] - 


8oy2bZ £ [v[ x - (1 - - i 2 !) (=0,0 + 
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■where m 1 is the tangent of angle of sweep of the midchord line, or 
<3M 

tan a 


4pV 2 s^jt 2 c r | 

?1 2 m*s 

1b + 5 sb; + 

1 - 

6* y 

|®0,0 

+ (i + 

4 m*s 
15 «c r 

(l | 16 m's , 4(1 - -X) 
\<o4 103 *c r 105* 

) a o^ 

+ (s 

1 m 1 s N 
3 Jtc r ^ 

) a l,0 + 

/ 1 2 m' £ 

j\ , /I 

, 8 

m's\ 


L - A 

(32 + 15 5tc 3 

:j%2 + 

105 

itc r ) B 'l,k ^ 

12* 

fe.jJl _ l} 

i*. + fen. 


. _aJ 1 

la , 


y 30* 64y 

' 2,2 \ 105* 

128j 

1 2,4 



i.-jjL + 

^5T] & o,2 


(C12) 


but 


H " 16 2 py2s ^ 6a 0,0 + ^1,0 + 4a 0,2 + ^l^ + 2a 0,4 + a l,^ 

p^sV 


8 


■( l6a 0,0 + ^1,0 + 4a 0,2 + 2a l,2 + 2a 0,4 + a l,4) 


The longitudinal center of pressure is at 

*cp\ = <2M _ 3 + Vo.g + *5*0,4 + *4 a l,0 + *5*1,2 + *6*1,4 + *7*2,0 + *8*2,2 + *9*2,4- 

jyfb r l6* 0 , 0 + 8^^ + l,a 0j2 + 2 *^ + 2a^ + 

(C13) 


where 



2 ale. + 

3 «c r 


1 A 

6* 


and kj), kj, ki,., k^, kg, ky, kg, and k^ are the coefficients 
of a Q ^ 2 , a Q ^, a 1>0 , a^ 2 , a^, a^ 0 , a^ 2 , and a^, 

in equation (C12). 


respectively. 
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TABLE I 

LOCATION AND MAGNITUDE OF WEIGHTING FACTORS TO 
REPRESENT CHORDWISE LOADING 


ChorcLwlse location, 
percent chord 

° 0 ,v 

Ol,v 


12.5 

0.2734 

0.0488 

0.0732 

37-5 

.1172 

.0762 

.0381 

62.5 ' 

.0703 

.0762 

-.0381 

87.5 

.0391 

.0488 

-.0732 














Wing details 


1.00 

U 

0.60 


Area, sq ft 
Aspect ratio 
Taper ratio 



Figure 2.- Details of wing-body configuration for wing A. All dimensions 

are in inches. 






1.857 



Figure 4.- Details of wing-tody configuration for wing B. All dimensions 

are in inches. 
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(a) M 0 a 0. 

Figure 5 .- Calculated sparnd.se loading. Wing Aj a - 4.0 1 










































(a) Twist = 0°. 

Figure il.- Comparison of calculated and experimental wing chordwise cen 

ter of pressure. Wing A. 





(b) Twist = ^.0°. 
Figure XI.- Concluded 







Figure 12.- Comparison of calculated and experimental wing chordwise 

center of pressure. Wing B. 
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Figure 13 - - 84- vortex pattern for wing showing location of nine control 

points . 
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